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1. Introduction 

Given the enormous success of phenomeno logical density dependent forces/functionals 
like the different variants of Skyrme, Gogny or relativistic mean fields in describing 
nuclear properties with mean field wave function [I], it was expected that going beyond 
mean field both restoring broken symmetries or considering configuration mixing would 
pave the way for a more satisfactory understanding of nuclear structure both in the 
qualitative and quantitative sides. Early attempts with the Skyrme functional [2] 
applying the Generator Coordinate Method (GCM) produced very encouraging results. 
Soon after, fully microscopic parity projection was implemented with the Gogny force 
[3] for octupole deformed states in a consistent framework and the comparison with 
other results as well as with experiment was encouraging. After that, many GCM and 
symmetry restoration calculations both with Skyrme and Gogny were reported. Here, I 
will only mention a few like the restoration of particle number symmetry that was tackled 
in [1] with the Gogny force. In this reference it was realized that a consistent inclusion of 
the direct, exchange and pairing terms of the interaction was needed for a divergent free 
evaluation of the energy. This was a real problem for functionals of the Skyrme type were 
the exchange and pairing fields of different terms of the functional were usually neglected 
(probably the most disappointing situation was the one of the Coulomb potential, as 
the consideration of its exchange and pairing fields increased by orders of magnitude 
the computational cost of the Skyrme type calculations). Recently, a regularization 
procedure to solve this problem and also to address the inherent self-energy and self- 
pairing problem of any energy density functional was proposed j5] and applied to the 
case of particle number projection (PNP) in [6]. Although the regularization procedure 
of [5] is basis-dependent and the impact of this arbitrariness has still to be checked, the 
results obtained so far [6] are very promising. 

Another problem that plagues beyond mean field calculations has to do with the 
prescription to be used for the density dependent part of the "Hamiltonian" / functional 
(quotes will be used through the paper to stress the dependency of the "Hamiltonian" 
operator with the intrinsic wave function through the density dependent part of the 
interaction). The first prescription used [2] was inspired by the three body force origin 
of some Skyrme functionals and is referred to as the "mixed density" prescription (see 
below). Although attempts to use other prescriptions are not uncommon El El CEO] 
it was assumed that the consistency requirements summarized in [TT] favored the use 
of the "mixed density" prescription. However, as pointed out in [121 EE] the "mixed 
density" prescription has a serious and in general insurmountable drawback when the 
density dependent part of the interaction depended upon non-integer powers of the 
density. As the mixed density is in general a complex quantity, non integer powers of it 
are not uniquely defined in the complex plane unless a "branch cut" is introduced and 
the Riemann sheet where the root lives is specified. In many cases, the mixed density 
is a real quantity and such a problem can be circumvented, but in general (including 
the most complex symmetry restoration of triaxial deformed states) it is complex and 
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the lack of a physically guided rule to choose the roots yields to projected energies 
that depend on external artifacts like the integration contours entering the definition 
of the projector operators p2J HB]- The only two known ways out to this problem 
in the general case are either to consider forces/functionals depending upon integer 
powers of the density (see below for a recently proposed functional) or to consider other 
prescriptions producing real densities. In the latter case, the most promising candidate is 
the "projected density" prescription [8], [91 [10] that, in addition, has some other appealing 
properties to be discussed below. In this paper, I am going to show by means of an 
example concerning parity projection that the "projected density" prescription yields to 
unphysical energies that favor infinitely deformed intrinsic states with diverging energies. 
This catastrophic behavior can also be extended to other spatial symmetry restorations 
like angular momentum projection or translational invariance. 

I would like to emphasize that the present discussion concerns the inadequacy of 
using on the same footing symmetrized densities and intrinsic ones to obtain the energy 
associated to the strongly repulsive density dependent interaction, which is characteristic 
of the nuclear Energy Density Functional (EDF). The present discussion has little to 
do with the solution [H] in terms of symmetrized density matrices of the symmetry 
dilemma of the Kohn-Sham Density Functional Theory (DFT). The symmetry dilemma 
refers to the use of KS-DFT density functionals which break in some circumstances the 
symmetries of the interaction and lead to Kohn-Sham wave functions that, in contrast 
to the full interacting wave function, cannot be assigned to an irreducible representation 
of the symmetry group of the interaction. 

The plan of the paper is as follows: in the next section, the "mixed density" 
prescription will be briefly discussed as a way to justify the introduction of the projected 
density prescription to be discussed in Section 3. In Section 4 I will try to convey what 
the open questions are regarding this issue as well as some thoughts about possible ways 
out. Conclusions will be given at the end of the paper. 

2. Mixed density prescription 

In the evaluation of the matrix element of a density dependent "Hamiltonian" between 
different states \<p) and \(p') a prescription for the density dependent part is needed. A 
common choice, proposed in [2] is the "mixed density" prescription that uses the overlap 



as the density of the density dependent term of the "Hamiltonian". In this definition 
p(f) = J2i=i S(r — fi) is the density operator. This prescription is based on the 
expressions obtained with three body forces in spin saturated systems. The use of 
this prescription makes the density dependent "Hamiltonian" to break the symmetries 
broken by the density and also makes it a complex (non-hermitian) quantity. One might 
think, after considering these two properties, that the mixed density prescription should 
be sentenced to oblivion. However, as was discussed in detail in Ref. [15], the energy, 




(1) 
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which is the physical quantity to be computed out of the above mentioned matrix 
elements, is a well defined quantity, i.e. it is invariant under the broken symmetries 
and it is a real quantity. To understand these facts let us consider briefly the case of 
rotational invariance. Let us now consider an intrinsic wave function \ip) that represents 
a deformed state and therefore is not characterized by the angular momentum quantum 
numbers. For standard (density independent) Hamiltonians the energy is independent 
of the orientation 

(cp\R+(n)HR(Q)\<p) ^ (<p\H\<p) (2) 

as a consequence of the scalar nature of the Hamiltonian, namely R(Q)H = HR(Q) for 
all Euler angles fl For density dependent "Hamiltonians" Hlp^)] it was shown in Ref 
[15] what the scalar property for overlaps below holds 

(tp\k{n)Hy>W) = (<p\H\p R }R(n)\<p') (3) 

where the shorthand notation pn = (ip\pR(£i)\(p') / (ip\R(ft)\ip') and up = 
(<p\R(Q)p\(p') / ((p\R(Q)\ip') is used. This property allows to extend Eq. (j2J) to the case 
of density dependent "Hamiltonians" showing that the ensuing energy is orientation 
independent. 

Concerning the reality of the energy, one has to realize that whenever a term 
like (iplHltpPtpiWip') appears in the expression of the energy, the symmetric term 
(tp'lHlpippWip) will also be present and weighted with the complex conjugate weight. 
It is not difficult to show [15] that the last term is the complex conjugate of the first 
and therefore the contribution of the two add up to a real quantity. This analysis 
is kind of a naive one as the density dependence in the Hamiltonian usually comes 
through a non integer power of the (complex) density what forces the introduction 
of "branch cuts" in the complex plane and to specify in which Riemann sheet the 
roots are defined. As a consequence of this difficulty it was shown in Ref, [13] (see 
also [12] for an early consideration of this issue) that the energy was not independent 
of the specific way it was computed. To be more specific, in the case of PNP with 
density dependent "Hamiltonians" consisting of a density dependent term with a density 
raised to a non integer power, the projected energy was dependent upon the integration 
path chosen in the complex plane to perform the integrations needed in the particle 
number procedure. This is a real challenge to existing theoretical models as most of the 
Gogny, Skyrme or relativistic models contain non-integer powers of the density in their 
definition. Recently, an exception has shown up in the market in the form of the so called 
Barcelona-Catania-Paris (BCP) functionals [161 E] which are only defined in terms of 
integer powers of the density and are therefore promising candidates to incorporate 
the self-energy renormalization techniques of [5] to produce well defined quantities in 
calculations beyond the mean field. 

A possibility that also comes up to solve the complex mixed density dilemma is to 
use the so called "projected density" prescriptions discussed next. 
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3. Projected density prescription 

The projected density prescription was postulated for Particle Number Projection 
(PNP) (HI |9] and Angular Momentum Projection (AMP) [10J to solve, in one hand 
, the problems associated with the complex nature of the mixed density prescription 
and, in the other, to give back to the density dependent part of the "Hamiltonian" 
the nice property of preserving fundamental symmetries like rotational invariance. In 
this way, the density dependent part of the "Hamiltonian" could be looked at as an 
effective part of a real Hamiltonian and not as a device to generate a non trivial density 
dependence in the energy functional [TO] . 

In the case of PNP the particle number projected density 

was considered in [H |9] and it was even proved [9] that this prescription as well as 
the mixed density prescription produced the same expressions for the Lipkin Nogami 
parameter in the strong deformation limit. This has been the prescription used in all 
particle number projected calculations with the Gogny force. In the rotational symmetry 
case the rotational invariant density corresponding to the J = projected state pj_g(r) 
was postulated [10] to be used in the density dependent part of the Hamiltonian. 

The "projected density" prescription looked reasonable and the only argument 
found against it in the PNP case [UJ was in the context of the GCM where, under 
some assumptions, it is possible to derive the Random Phase Approximation (RPA) as 
a small amplitude approximation of the GCM [18]. This idea can be extended to the 
case of density dependent interactions [UJ, [19] with the "mixed density" prescription to 
produce the RPA with the same matrix elements that appear in the standard derivation 
of the RPA for density dependent forces [20]. However, this is not the case for the 
"projected density" prescription or other prescriptions postulated in the literature [7]. 
It has also to be mentioned that the "projected density" prescription produces in the 
context of the PNP and in the strong deformation limit an approximate expression for 
the projected energy (Kamlah expansion, see [2T] for a detailed discussion) including 
a parameter, to be identified later with the Lagrange multiplier of the particle number 
constraint, that includes the rearrangement term characteristic of mean field theories 
with density dependent forces [H [9], [TTJ. This property do not hold, however, for 
spatial symmetries like rotations as the leading term of the Kamlah expansion for the 
projected density is not the intrinsic density as in the case of PNP but some kind of 
symmetry preserving average (see below for the parity symmetry). Therefore, in the 
general case the density dependent interaction of the underlying mean field is different 
from the density dependent interaction to be used in projection (i.e. both densities 
are "decoupled" in a sense to make clear below) suggesting the possibility to use the 
"projected density" prescription also at the mean field level. That would make the 
Hartree-Fock-Bogoliubov (HFB) equations more involved due to the new rearrangement 
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Figure 1. (Color online) On the left hand side various energies for the nucleus 36 Ar 
at a fixed value of fi% of 0.23 are plotted as a function of the octupole moment Q3. The 
blue dashed increasing curve corresponds to the intrinsic HFB energy. The red full line 
is for the projected energy to positive parity using the "mixed density" prescription. 
The dotted green down-sloping curve is for a calculation of the HFB energy where 
the intrinsic density of the density dependent interaction is replaced by the projected 
density to positive parity. Finally, the purple full line down-sloping curve is for the 
projected energy corresponding to positive parity and with the parity projected density 
in the density dependent part of the interaction. The result was obtained with the 
Gogny D1S force. On the right hand side, I plot the same quantities as before but for 
the heavy actinide and octupole deformed nucleus 224 Ra. 



term but no fundamental problems were foreseen. As I will show below this possibility 
is also doomed to fail. 

In order to show the catastrophic consequences of having in the density dependent 
part of the "Hamiltonian" a quantity that is decoupled from the intrinsic density I will 
consider the case of parity projection of octupole deformed intrinsic states. The Gogny 
[22] force D1S [23] will be used for the calculations and I will proceed along the same 
lines as our early study of parity projection with the Gogny force of Refs. (3j El]. The 
parity projection operator to parity n = ±1 is given by 

P* = ^={I + *fi) (5) 

where fl is the parity operator. The projected energy is given by 

E _ (<j>\HP*\<j>) _ (^gM+7T^g%) 

n (0|P^|0) " (0|0) + vr(0|fl|0) 
The definite expression for the projected energy depends upon the prescription used for 
the density dependent part of the "Hamiltonian" . For the "mixed density" prescription 
the standard intrinsic density p(r) has to be used in the density dependent part of the 
Hamiltonian for the evaluation of its mean value (first term in the numerator of Eq. (JH]) 
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) and the density 9{r) = ^jgrP in the evaluation of (0|M|0)/(0|n|0). In the case 
of the "projected density" prescription the density 

. =+ i (01^/3(0^10) p(r) + vr(0|n|0)g(r) 

PpRQjU (0IP+4I0) (l + 7r<0|%>) 

will be used in the density dependent part of the "Hamiltonian" . In this expression 
the averaged matrices p{f) = | (p(f) + p(— r)) and #(f) = | (6*(f) + 0(— r)) have been 
introduced. In the strong deformation limit (corresponding to (0|II|0) — > 0) the 
projected density goes to the average density p(f) instead of the intrinsic one. This 
is the fact that is at the heart of the problems encountered below for the projected 
energy. 

The positive parity projected energy E + has been computed using the above 
formulas and prescriptions for given intrinsic configurations with different octupole 
moments and in different nuclei. A light nucleus 36 Ar has been chosen as an example 
of a system with a density which is mainly surface and also as a system where the 
intrinsic ground state is reflection symmetric (i.e. positive parity in this case). On the 
other hand, the heavy nucleus 224 Ra has been chosen as an example of an octupole 
deformed nucleus and with a density which is mostly volume. The results obtained 
are shown in Fig. [1] as a function of the octupole moment Q%. Several energies are 
computed: the mean field energy computed in the standard way E, the mean field energy 
computed with the parity projected density in the density dependent "Hamiltonian" 
Epjy (the modified HFB framework discussed above), the positive parity projected 
energy computed with the "mixed density" prescription E an d finally the positive 
parity projected energy computed with the "projected density" prescription E + pj> 
First, it is observed the typical parabolic behavior of E as a function of Q3 for a non 
octupole deformed nucleus like 36 Ar as well as the presence of a minimum for the case 
of the octupole deformed 224 Ra case. The positive parity projected energy E.y^ 
obtained with the "mixed density" prescription has a behavior close to the mean field 
energy being lower than it in all the cases (see [3]) and quickly converging to the mean 
field energy for increasing values of the octupole moment. This is a consequence of 
the behavior of the mean value of the parity that quickly decreases to zero as the 
octupole moment increases reaching the strong deformation limit. On the other hand, 
the two corresponding energies computed with the "projected density" prescription show 
a very pronounced decrease of tens of MeV that keep going as the octupole moment is 
increased. In these two cases, the configuration minimizing the corresponding energies 
would be an infinitely deformed octupole configuration with energies going to minus 
infinity. Obviously this is a catastrophic situation that will lead the system to unphysical 
energies and intrinsic configurations. Therefore the "projected density" prescription, at 
least for spatial symmetries (see below), can be categorized as unphysical and its use has 
to be avoided. Concerning the particle number projected density prescription widely 
used with the Gogny force, the reasonable results obtained so far clearly demonstrate 
that no catastrophe is taking place in that case. However, and from a fundamental (as 
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Figure 2. On the left hand side, the contour plots of the intrinsic density of the 
nucleus 224 Ra with an octupole moment of 4 6 3 / 2 (roughly corresponding to is ground 
state) are depicted. The contour lines start at p = 0.02fm~ 3 for the outer contour and 
end at p — 0.14fm -3 for the inner one. On the right hand side, contour plots of the 
positive parity projected density corresponding to the intrinsic wave function of the 
left hand side plot. The contours are the same as in the left hand side plot. 



opposed to phenomenological) point of view, the use of a prescription which yields to 
catastrophic results when used in other contexts should also be avoided. 

In order to understand the catastrophic results obtained it is convenient to analyze 
the structure of the density dependent energy which is given in the case of the Gogny 
force by the following integral 

(<p\Hdd\<p) = | £(! + y - (*o + ~)<V) / d 3 fp T p T ,p a (8) 

t,t' ~~ ~ 

where p T are the intrinsic densities with isospin r. The replacement of the standard 
p a {f) by Pp ro j(^*) produces a "mismatch" in the integrand because the new density 
pp ro j does not have the same shape as p as can be seen in Fig. [2] where contour plots of 
the two densities are plotted for the nucleus 224 Ra at an octupole moment Q 3 = 4b 3 ^ 2 
(which corresponds to the ground state of this octupole deformed nucleus). As pp ro j 
is reflection invariant it must have "humps" at each side of the nucleus whereas the 
intrinsic density only has a "hump" at one of the sides (which one depends on the sign 
of the octupole moment). The two "humps" of the projected density will have half 
the density of the one of the intrinsic density as the projected density is roughly given 
(strong deformation limit) by the average density p(r) = \ (p(r) + p(— r)) (see Eq. ([6])). 
As a consequence of the mismatch of the different factors, the integrand corresponding 
to the "projected density" prescription is slightly reduced as compared to the integrand 
with the intrinsic density in the density dependent term. As a consequence, the value of 
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the integral is also reduced producing a reduction of the (repulsive) density dependent 
energy that is magnified by the large values of the density dependent strength parameter 
t 3 ( 1562.22 MeV fm 4 , 1609.46 MeV fm 4 and 1350.0 MeV fm 4 for Gogny DIM [26], 
D1N[25J and D1S [23j, respectively). 

The mismatch between the projected density and the intrinsic density will also be 
very pronounced in the case of quadrupole deformation and associated to rotations 
in space as well as in the case of spatial translations. It is therefore evident that 
those two cases will also be plagued with the catastrophic decrease of the projected 
energy for infinitely deformed configurations. Therefore in the case of angular 
momentum projection and restoration of translational invariance the "projected density" 
prescription is doomed to badly fail. 

4. Open problems 

From the results presented here it is clear that, at least for spatial symmetries, the 
only possibility to get rid of all the different problems encountered in the extension of 
the density dependent forces paradigm to calculations going beyond mean field is to 
solely consider forces/functionals depending upon integer powers of the density. Apart 
from some scarce versions of Skyrme, there is a recently proposed functional with this 
property, the Barcelona-Catania- Paris (BCP) functional [TU [17] (see also an article in 
this issue by Baldo et al concerning this functional [27J) which is based on a fit to a 
realistic nuclear matter equation of state followed by a Local Density Approximation 
map to produce a functional for finite nuclei. Additional finite range surface terms, 
spin-orbit, Coulomb and density dependent pairing are ingredients of this functional. 
Given the good performance of BCP at the mean field level [TBI [T~7l 127] . this functional 
can become a serious contender for beyond mean field studies in the spirit of the Energy 
Density Functional (EDF). 

Another alternative to get rid of the above problems would be to use approaches 
based on the Gaussian Overlap Approximation (GOA) [21] or its more complete 
accurate and satisfactory version (the Topological GOA, see for instance [28] for a 
recent treatment) to perform the symmetry restorations and/or configuration mixing 
calculations (in the GCM spirit). This possibility has been explicitly suggested in [29] 
in the framework of the Bohr Hamiltonian method but the philosophy behind was also 
implicit in early applications of the method [30]. The reason why this proposal will 
cure the problems discussed before is because the residual interaction is defined as 
second derivatives of the energy (where the density dependent term is well defined) and 
therefore are free from the above problems of divergences and complex valued quantities. 
Obviously, these ideas can be generalized to a framework where overlaps are computed 
exactly but the residual interaction is defined as a second derivative in a way that closely 
resembles the RPA for density dependent forces/functionals. Work along this direction 
is in progress [T9] . 
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Conclusions 

I have explicitly shown with an example using parity projection that the "projected 
density" prescription, advocated by some authors to overcome some difficulties 
encountered in the application of the "mixed density" prescription in applications 
beyond the mean field, leads to catastrophic consequences for the projected energy as 
this prescription favors a system where the reflection symmetry violation is maximal and 
the projected energy in unbounded by below. This consideration can also be extended to 
other symmetry restorations involving spatial symmetry transformations like rotations 
and translations. The result is the consequence of the mismatch between the shape 
of the intrinsic and the projected densities and the strongly repulsive character of the 
density dependent term. 
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